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Abstract: The equation of state of a family of alternative linear chiral models in the mean 
. field approximation is discussed. We investigate the analogy between some of these models 

' with current models in the literature, and we show that it is possible to reproduce very well 

, the saturation properties of nuclear matter. 

OO 

^ ' Since long time ago it has been known that chiral symmetry (CS) is essential for the 

understanding of hadronic reactions. On the other hand its role in nuclear physics was 
CJ I much less clear. The simplest implementation of chiral symmetry in this context is the 

^ . linear sigma model 1^, where the scalar meson plays a dual role as the chiral partner 

of the pion and the mediator of the intermediate-range nucleon-nucleon attraction. 
However, this model was not able to reproduce the saturation properties of nuclear 
matter at the mean field level and, therefore, was left in oblivion and much attention 
was given to field theoretical models like QHD in its inumerous variants. During the 
eighties some important works tried with relative success to cure the problems of the 
linear sigma model through the introduction of a vector-scalar coupling [^j. However, 
these improved versions of the model were shown to fail in giving a good description 
of finite nuclei properties. The hope, and to some extent increasing evidence ||^, that 
chiral symmetry is really important for nuclear physics motivated the implementation 
of nonlinear realizations of CS, some of them being very successfuU in the applications 
to finite nuclei p, ^. 

Alternatively, with the introduction of a scalar glueball field, a reasonable descrip- 
tion of nuclear matter and closed shell nuclei was achieved in the context of the linear 
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sigma model [Q . The idea behind the introduction of the glueball field was to create an 
effective lagrangian that respects the low-energy theorems of scale and chiral symme- 
try. The authors of ref. [§] have also noticed that the coupling of the vector field with 
the glueball field rather than the sigma field is essential to provide a good description 
of finite nuclei, in agreement with the conclusions in . 

Prospectively, in this paper we show that it is possible to get a good description 
of nuclear matter and eventually of finite nuclei with a generalized formulation of the 
linear sigma model, that contains only the original scalar, pseudoscalar and vector 
meson fields without additional couplings between the mesons. To generate the model 
we use, as motivation, the quark meson coupling model (QMC) which is a model 
for nuclear matter and finite nuclei in which the quark structure of the nucleons is 
explicitly considered. The basic assumptions of the model are that effective meson 
fields couple directly to quarks confined in non-overlaping bags, and that nucleons obey 
the Dirac equation in the effective meson mean fields whose sources are the quarks in 
the bag. The QMC model shares important features with the Walecka model (QHD), 
such as the relativistic saturation mechanism. The final equation for the total energy 
per nucleon in nuclear matter is identical to that of QHD. The effect of internal quark 
structure enters only through the effective mass and the self-consistency condition for 
the a field. The self-consistency condition, on the other hand, differs from that of QHD 
through a scalar density factor. This is equivalent of having a coupling constant which 
depends on the scalar field and consequently on the density. A similar dependence is 



also found in the chiral confining model of ref. |10]. In that work it is also shown 
via a Dirac-Brueckner-Hartree-Fock calculation that significant changes in the nuclear 
matter saturation may be observed when the mesonic coupling constants depend on 
the density. Another kind of density dependent coupling constant was demonstrated 



to come from relativistic SU (6) model |TT 



Since it is well known that models with linear realization of chiral symmetry have 
difficulty to saturate the nuclear matter we pose the question whether a density depen- 
dence of the mesonic coupling constants, claimed in the above mentioned works, could 
provide a reasonable saturation mechanism in this class of chiral models. Therefore, 
we will allow the coupling constants in the linear sigma model to be dependent of the 
scalar field. We start with the following lagrangian: 
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-\f,^F^^ + \mlV,V^ - \x{a' + tt^ - u^^f , (1) 

where F^i, — — 9,^V^. Here ^, V^, a and tt are respectively the nucleon, vector 
meson, scalar meson and pion field, u is the vacuum expectation value of the sigma 
field which gives mass to the nucleon. In this version of the sigma model there is no 
symmetry breaking term and the symmetry is realized in the Nambu-Goldstone mode. 
The effective couphng constants are given by 

' (l + ^^)" ' 



where is the pion- nucleon coupling contant, which value is fixed at the experimen- 
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tal value: |^ = 13.6, and g^, the vector meson- nucleon coupling constant, is a free 
parameter, a and /3 are constants which define different families of models. After 
the introduction of field dependent couplings the final functional dependence of the 
lagrangian density on the fields is very different from the original linear a model. How- 
ever we preserve the linear realization of the chiral symmetry. This is explicit in the 
new factors which depend only on the combination o"^ + tt^. 

We now shift the scalar field defining a = u + s , such that (s) = in the vaccum. 
The lagrangian becomes: 

C = n^^izd^ - glVn -M*- igl^.f ■ 7f]vl/ + ]^{d^sd^'s - my) 

+ -dj ■ 5^7f - -F.,F^- + -mlV.V - ^s{s^ + tt^) - ^{s^ + tt^ (4) 
2'^ 2 2 2M ^ ' 8M2 ^ ' ^ ' 

where = 2\v? ^ M — g^^u and M* — ^{M + g-^s). The non vanishing vaccum 
expectation value of the a field generates both the cr and nucleon masses. In a mean 
field (MF) approximation the energy density is given by 



(5) 



where = , = , is the barion density: 

kp 

J d'k , (6) 
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and 7 is the spin-isospin degeneracy (7 = 4 for nuclear matter). 

To write Eq. we have used the fact that in MF we can write the effective couphng 
constants as 
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The effective nucleon mass, M*, is determined by minimizing Eq. (^) with respect 
to M*: 

4/3+2-4a 
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At this point it is interesting to compare the MF resuh of our model with models 
existing in the literature. In particular, in ref.|T2| a unified model Lagrangian density 
was proposed to connect QHD with the Zimanyi-Moszkowski (ZM) model |13| and a 
modified version of the ZM model, called ZM3 in ref. |]T^. In the MF approach, the 
energy density of this unified model is given by 

kp 
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where the gap equation for the effective nucleon mass is 



2Cl \ {M*/Mf 
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The values of a and h for the different models given in Eq. ([l0|) are given by 



model 


a 


b 


QHD 








ZM 





1 


ZM3 
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1 



Comparing Eqs.(^ and (0) and Eqs.(^ and (|TT]) we can see that in MF our model 
can reproduce QHD, ZM and ZM3 models if a and /3 assume the values given bellow. 



provided that we redefine our as C^. 
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model 


OL 


p 


QHD 


-1/2 





ZM 


3/2 





ZM3 


3/2 


1 


linear-o" model 









We have included in the table above the values of a and /? that reproduce the original 
linear-o" model [|l], H, |15|. As it is well known, the original linear-cr model does not 
reproduce the nuclear matter ground-state properties. This happens because in this 
model M* goes to zero at small values of the nuclear density and this generates a very 
strong attraction. In our model this excess of attraction is avoided by the decrease of 
the scalar coupling constant as a function of the density, as can be seen in Eq.(|^) (since 
M* is a decreasing function of the density) for a < and a > 1/2. 
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Figure 1: The scalar meson mass as a function of a for different values of (3. The circle and the 
cross stand for the ZM and ZM3 results respectively. 



We will analyze our model in the region a > 1/2 and for three choices of (3: (3 = 0, 
P = I and P = a. The values of and are fixed to reproduce the energy per 
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nucleon E/N — —15.75 MeV at kp — 1-3 fm~^. Since is fixed at its experimental 
value, different values of are related with different values of m^. In Fig.l we show 
the value of the scalar meson mass obtained from that saturates the nuclear matter, 
as a function of a. The particular results obtained in the ZM and ZM3 models are 
assigned in this figure by a circle and a cross respectively. Prom this figure we can 
see that the choice j3 = a provides very smaller values for and delimits the region 
where the symmetric scaling {a = j3) is valid. If one thinks about finite systems, 
is associated with the range of the scalar interaction. Therefore, smaller values of 
provide longer ranges of the interaction in the configuration space. 

900 
800 
700 
600 

> 

^ 500 

* 

^ 400 
300 
200 
100 

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 

a 

Figure 2: The effective nucleon mass at saturation density as a function of a for different values of 
/3. The circle and the cross stand for the ZM and ZM3 results respectively. 

In Figs. 2 and 3 we show the effective nucleon mass at saturation density, Mq, and 
the incompressibility, as a function of a. As we can see by these figures, a large 
region of values for K and Mq can be generated by these models and we can choose the 
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values of a and /3 that can better reproduce the phenomenology. It is also interesting to 
notice that the values generated for the incompressibility by these models are generally 
small, as compared with QHD for instance. 

There is a strong correlation in relativistic models between the effective nucleon 
mass at saturation density, Mq, and the strenght of the spin orbit force in nuclei |T6| , 
PT| , |TB[. Typical values of Mq/M in successful models are roughly 0.6, which yield spin- 
orbit splittings close to experimental values. Accurate reproductions of observables in 
finite nuclei, including the spin-orbit splitting in particular, tightly constrain the value 
of Mq/M in the range 0.58 < Mq/M < 0.70 [|T^. Since in linear chiral models in 
general Mq/M is too close to unity, this was found as one of the reasons for the 
failure of the models considered in ref.[Q to give a good description of finite nuclei. As 
can be seen by figure 2, allowing the coupling constants to decrease with the density 
we can generate a linear chiral model which gives Mq/M in good agreement with 
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phenomenology. We can also see that the choice /3 = does not go in this direction 
generating M^/M very close to unity, indicating its almost non-relativistic character. 
The choice a = P and P = I (for a around 1.5 (ZM3)) seems to be a good candidate. 
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Figure 4: — S" as a function of a for different values of (3. 



Another qualitative ingredient associated with the spin-orbit splitting in finite nuclei 
is the difference between the depth of the vector and scalar potencials, V — S {S = 
M* — M, V = g*'^pB/ml). The higher the difference the higher the splitting. In Fig. 
4 we present the difference V — S as a function of a. Typical values of this quantity, 
that give reasonable values for the spin-orbit splitting, lies usually in the region from 
550 to 650 MeV In our models this region can be obtained for the cases a = (3 and 
/3 = 1 in agreement with the analysis done above. Of course this is only an indication 
that good results can be obtained for finite nuclei. A calculation for finite nuclei with 



these models is under way |20 
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To summarize, we have proposed a new class of models for nuclear systems where 
chiral symmetry is implemented in a linear realization. No new degress of freedom or 
new interactions other than that present in the original linear-a model |2T | 



were 



needed. In these models loops are not an essential ingredient to obtain saturation of 
nuclear matter as in refs. pT| , p!5[] . We found that this class of models reproduces very 
well the expected saturation properties of nuclear matter. In a mean field approxima- 
tion and for some particular values of the parameters, these models can reproduce the 
well known derivative coupling ZM models and QHD. Therefore ZM models and QHD 
can also be understood as reminiscents of linear chiral models with couplings given 
by Eqs. (0) and @, after a mean field approximation and proper choice of a and /3. 
Let's remark that ZM models have a desirable nonrelativistic limit in the saturation 
mechanism. We are aware that for finite nuclei, surface terms can not be neglected 
and the results for nuclei spectra may differ from the original models |jl2|. Finally we 
would like to point out that to the best of our knowledg it is the first time that this 
class of linear chiral models can reproduce well nuclear matter phenomenology. 
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